Magnetoconductance of the Corbino disk in graphene 
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Electron transport through the Corbino disk in graphene is studied in the presence of uniform 
magnetic fields. At the Dirac point, we observe conductance oscillations with the flux piercing the 
disk area characterized by the period $o = 2 (h/e) ln{Ro/Ri), where Rc, (Ri) is the outer (inner) 
disk radius. The oscillations magnitude increase with the radii ratio and exceed 10% of the average 
conductance for Ro/Ri ^ 5 in the case of the normal Corbino setup, or for Ro/Ri ^ 2.2 in the case 
of the Andreev-Corbino setup. At a finite but weak doping, the oscillations still appear in a limited 
range of ^ 'E'™^") away from which the conductance is strongly suppressed. At large dopings 
and weak fields we identify the crossover to a normal ballistic transport regime. 



An atomically thin carbon monolayer (graphene) is 
widely considered as a successor of silicon in future elec- 
tronic devices [T] . Investigations of the low-energy prop- 
erties of graphene, governed by the massless Dirac equa- 
tion, constitute new and thriving sub-area of condensed 
matter research |2] . Particularly striking feature of clean, 
undoped graphene samples is that zero density of states 
is accompanied by a nonzero, universal value of the con- 
ductivity 4e^/(7r/i) [3H6]. This is a basic signature of 
the so-called pseudodiffusive regime, in which transport 
properties of graphene are indistinguishable from those 
of a classical diffusive conductor [?]• In this regime, the 
applied magnetic field does not affect the conductivity 
[SI H] and higher current cumulants [TU]. Prada et al. 
also show that for high dopings and magnetic fields, the 
pseudodiffusive behavior is recovered at resonance with 
the Landau levels (LLs) in the absence of disorder. 

Numerous studies of graphene magnetoconductance fo- 
cus on nanoribbons [11] , Aharonov-Bohm rings [12l [l3] , 
antidot lattices [14], and weak- localization effects in 
chaotic nanosystems [T31 [TB] . Cheianov and Fal'ko [T7] 
showed the conductance of a circular p-n interface is in- 
sensitive to the weak applied field. The author, Recher, 
and Wimmer recently identified the crossover from the 
pseudo-diffusive to the quantum-tunneling regime |18j . 
which is characterized by a power-law decay of the con- 
ductance G (X L~°' (where L is the length of a sample 
area and a is a geometry-dependent exponent) and ap- 
pears for quantum billiard in undoped graphene at zero 
field. In the case of the Corbino disk with the outer ra- 
dius i?o and the inner radius Ri (see Fig. [T]) we have 
L = — Ri and a = 1, leading to the reciprocal de- 
cay of G for Ro ^ Ri- As the tunneling regime shows 
up generically for billiards having (at least) one narrow 
opening [18 , the discussion of magnetic field effects — at 
least on a basic example — is desirable. 

In this Rapid Communication, we analyze theoretically 
magnetoconductance of the Corbino disk in graphene at 
arbitrary dopings and magnetic fields. The paper is or- 
ganized as follows: We start from the mode-matching 
analysis for the disk attached to heavily-doped graphene 
leads, which employs the total angular momentum con- 
servation in a similar way as early works employed trans- 



verse momentum conservation for the strip geometry 
^ [4] . Then, we discuss separately the zero- and finite- 
doping situations, and present the system phase dia- 
gram in the field-doping parameter plane. The findings 
of Refs. [1^ for the pseudo-diffusive regime are repro- 
duced for Ro/Ri S 2. The novel feature is a periodic 
(approximately sinusoidal) magnetoconductance oscilla- 
tion, visible in undoped or weakly doped disks with larger 
radii ratios, and recovered at LLs for high dopings. Fi- 
nally, we extend the analysis to the normal-graphene- 
superconductor (Andreev-Corbino) setup. 

The analysis starts from the Dirac Hamiltonian in 
a single valley [TU] , which is given by 

H =:VF{p + eA)-(T + U{r), (1) 

where vp = lO^m/s is the Fermi velocity, cr — {ax,<Jy), 
p = —ih{dx, dy) is the in-plane momentum operator, 
the electron charge is — e, and we choose the symmet- 
ric gauge A — ^{~y,x). The electrostatic potential 
energy U{r) — Uq in the disk area {Ri < r < Ro), 
otherwise U{r) = Uoo- Since the Hamiltonian ([T]) com- 
mutes with the total angular momentum operator — 
—ihdip -t- haz/2, the energy eigenfunctions can be chosen 




FIG. 1: The Corbino magnetometer in graphene. The current 
is passed through the disk-shaped area with the inner radius 
Ri and the outer radius Ro in a perpendicular magnetic field 
B = (0,0,-6). The leads (yellow) are modeled as infinitely- 
doped graphene regions. The gate electrode (not shown) is 
placed underneath to tune the doping in the disk area. 
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as eigenstates of 

^,(r,^)=e^(^-V2).(^^X,TM^ (2) 

where j is an half-odd integer, s =t, i denotes the lat- 
tice pseudospin, and we have introduced the polar co- 
ordinates {r,ip). The Dirac equation now reduces to 
HjXjir) = Exj{r), where Xjir) = [Xjtir)'Xj-l{r)f , and 



-lhVF<Jxdr + U{r)^ 



eBr 
2h 







3 + 1/2 eBr 



(3) 



Subsequently, the scattering problem can be solved 
separately for each j-th angular momentum eigenstate 
incoming from the origin (r = 0). As the angular depen- 
dence of the full wavefunction ([2| does not play a role for 
the mode matching, the analysis limits effectively to the 
one-dimensional scattering problem for the spinor Xj (f) ■ 
We model heavily-doped graphene leads by taking the 
limit of U{r) = Uqo — >■ ^oo (hereinafter, the upper sign 
refers to the conduction band, and the lower sign refers 
to the valence band) , and define the reflection (transmis- 
sion) amplitudes Vj (tj). For the inner lead (r < i?i), the 
wavefunction can be written as 



-1 



(4) 



where the first term represents the incoming wave, and 
the second term represents the reflected wave. We further 
introduced kao = \E — Uoo\/{fi'Vp) oo. For the outer 
lead (r > R^) the wavefunction is 



(o) + 



(5) 



and represents the transmitted wave. Defining fco = \E — 
\jQ\lihvF) for the disk area (_Ri < r < i?o), we write the 
wavefunction in a similar form as considered by Recher 
et al. [20j for the eigenvalue problem, namely 



where z^^i = [2(j + Sj)]-^*^ , Zj^2 = 2{l3/k^y^+^/'^ (with 
Sj = ^sgnj, (3 = eB/{2h)), and 

^^""^ I t/(a,.,7,.,/?r2), ^ = 2, 

with \ for s =t,l, "js = \ W-s + \ls\ + 1) - 

/cq//?], and = |/s| + 1- M{a,b,z) and U{a,b,z) are 
the confluent hypergeometric functions [21j . Solving the 

matching conditions xj'^(i?i) = x^^^(i?i) and x'f\Ro) — 



xI'^\Ro) we find the transmission probability for j-th 
mode 



16(fcoV/5) 



|23-1| 



r(7: 



.r(ajt) 

where T{z) is the Euler Gamma function, and 



(8) 



w 



j(2) 



Without loss of generality, we choose B > 0. For B < 
one gets T^^B) = T^j{-B). 

First, we consider the zero doping limit, for which Eq. 
([8| simplifies to 



T,{ko^O) = 



1 



cosh^ [£(j + $d/$o)]' 



(10) 



where C = ln(i?o/i?i), $d = 7r(i?2 - i?f)B is the flux 
piercing the disk area, and $0 = 2 {h/e)£. We observe 
varying the ratio $d/<i>o affects Tj{ko — >■ 0) similarly as 
changing boundary conditions affects the corresponding 
formula for the strip geometry [3l |4] . (Notice that Eq. 
(10 1 is insensitive to the flux piercing the inner lead.) 



The disk conductance follows by summing over the modes 
G^goJ2T,{ko^O)=f2G,^cos('^^^), (11) 

where go — 4e^//i is the conductance quantum (the factor 
4 includes spin and valley degeneracy), and the Fourier 
amplitudes are 



Go 



4^^(-)"'mgo 
£^ sinh(7r2TO/£) 



(m>0). (12) 



The conductance given by Eq. (11) shows periodic oscil 



lations with the average value Go equal to the pseudo- 
diffusive disk conductance [18]. (Thus, the averaging over 
$d for the disk corresponds to the fictitious averaging 
over boundary conditions for the strip.) The approxi- 
mate formula G w Go -1- Gi cos(27r$d/$o) reproduces the 
full expression with the 1% accuracy for Ro/Ri ^ 10. 
The oscillations magnitude AG = Gmax — Gmin ~ 2|Gi| 
converges rapidly to with Ro/Ri — > 1 (the pseudo- 
diffusive transport regime), in agreement with earlier 
works [SHIO] reporting no field dependence of the con- 
ductance. For instance, we obtain AG < 4-10~^Go for 
Ro/Ri ^ 2. In the tunneling regime, the oscillations 
magnitude of AG > O.IGq is reached for moderate radii 
ratios Ro/Ri ^ 5. In this regime, when f^d/^a is half-odd 
integer, the major contribution to the conductance orig- 
inates from a single mode (T_$^/$g = 1) and we have 
G = Ginax (with Gniax 90 for i?o > ^i)- For other 
values of the conductance is generally dominated by 
the two modes, with j± = — int(<I>£;/$o "F |) — |, which 
became equivalent for ^d/^o integer, when G = Gmin 
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FIG. 2: Conductance as a function of the magnetic field and 
the doping for Ro/Ri = 10. (a) Magnetoconductance at weak 
doping (specified for each curve by koRi = 10"'' 10~\ with 
ko — \E — Uol/hvp)- The zero-doping magnetoconductance 
is also shown (red line), (b) Magnetoconductance at large 
doping (fco-Ri = 0.2 0.6). (c) Conductance as a function of 
doping at fixed magnetic field (specified by the flux piercing 
the disk area $d = 0~ 4$o, with $o = 2 (h/e) In Ro/Ri.) 



(and GnunRo/Ri ^ 
zero-field situation 



8go for Ro ^ Ri) reproducing the 



We now complement the discussion by analyzing a fi- 
nite doping case, to find out how stable are the con- 
ductance oscillations when the gate voltage is controlled 
with a finite precision. For fco > Eq. ^ is well de- 
fined for arbitrary j provided that |fcg//3 = (kalg)'^ /2 ^ 

n = 1,2,... (LLs) with Ib = \/h/eB the magnetic 
length. In such case, the asymptotic form for large fields 

"[-co- 



is T, 



(n) 



cosh 



2n + $d/$o) 



leading to con- 
ductance oscillations as obtained above, see Eq. (11). In 
fact, Tj(fco 0) given by Eq. (10 1 are reproduced for 
n = 0, showing the conductance oscillations for an un- 
doped disk can be rationalized in terms of resonant trans- 
port through the zero-th LL pinned at the Dirac point. 

The results for G, obtained by numerical summation 
of Tj-s given by Eq. (|8| for Ro/Ri — 10, are shown in Fig. 
[2j We first compare, in Fig. [2ja), the zero-doping mag- 
netoconductance (red line) given by Eq. (11) with those 
obtained for dopings varying from fco^i ~ 10"* to 0.1, 
with the steps of one order of magnitude. Weak-doping 
curves follow the zero-doping one for first few periods, 
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FIG. 3: Phase diagram representing the tunneling, field- 
suppressed and ballistic transport regimes in the field-doping 
parameter plane. Solid lines corresponds to G = Go for the 
radii ratio Ro/Ri = 10. Dashed lines depict borders of the 
tunneling ($d < 4'™'"', red line) and ballistic (2rc > Ro — Ri, 
blue line) transport regimes. (Notice the logarithmic scale for 
ko{Ro~Ri) < 1.) Inset shows the crossover into the tunneling 
behavior for the first Landau level (n = 1) in the magnified 
horizontal scale (5„ = ^kgl^ — n). 



when 



<^d<<fT''^—HkoRi). 

e 



(13) 



For higher fields, G decays as e~(^°~^')^/(^'^^. The high- 
G area, limited by Eq. ( 13 ) shrinks rapidly with increas- 
ing ko- However, for high dopings {ko{Ro — Ri) ^ tt) 
the results start to follow the semiclassical picture, sim- 
ilarly as for the two-dimensional electron gas (2DEG) 
P^l . The high-G area expands with ko (see Fig. [2]d), 
as it is now limited by the condition 2rc > Ro — Ri 
(with Tc = fco^g the cyclotronic radius), characteriz- 
ing the ballistic transport regime. In particular, for 
2rc > Ro + Ri we have G/go « 2koRi, approaching the 
result for zero field [HI [23] . At high magnetic fields (for 
which 2rc < i?o ~ ^i) we enter the field-suppressed trans- 
port regime, in which G ^ g"*^^""^'^ /(2'i,) again, except 
from the isolated peaks (see Fig. [2j^c) for the plot in a log- 
arithmic scale) , which correspond to the resonances with 
LLs, and shrink with the field in the absence of disorder. 
At each resonance, the zero-doping field-dependence of 
G is approached for the high field. 

The behaviors described above are presented in a con- 
densed form in the phase diagram shown in Fig. [Sj Col- 
ored areas represents the regions in the field-doping pa- 
rameter plane where G > Gq (with the borders G = Go 
marked by solid lines). We also show (with dashed lines) 
the limiting values of the magnetic field, at which the 
crossovers from the tunneling (left) and from the ballis- 
tic (right) to the field-suppressed transport regime occur. 
For the first LL, we demonstrate in a quantitative man- 
ner (see the inset) how, starting from the ballistic regime 
and enlarging the field (but keeping (A:o's)^/2 ~ 1), one 
restores the tunneling behavior, characterized by a chain 
of isolated islands of G > Gq on the diagram. 
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So far, we have considered the disk attached to normal- 
metaUic leads. For the Andreev-Corbino setup, with one 
normal and one superconducting lead, the conductance 
is expressed in terms of Tj-s given by Eq. Q as [23] 



G 



NS 



230 E 



rp2 



(2- 



(14) 



In particular, G^^ is still a periodic function of at the 
Dirac point, and its Fourier decomposition G^^($d) — 
+ cos(27rm$rf/$o) leads to 

G^^iC) = G„{C), Gl^iC) = 2G„,{2C). (15) 

Although we have G^/G^ 1 for Ro/Ri oo (and 
any to), magnetoconductance oscillations are noticeably 
amplified for moderate radii ratios. For instance, the 
magnitudes AG/Gq > 0.1 are now reached for Ro/Ri ^ 
2.2. At finite dopings, Eq. ^ for ^Jf^'^ holds true, and 
the phase diagram in the held-doping parameter plane 
(Fig. [3]) is almost unaffected. We further notice that for 
available ballistic graphene samples 2i?o < 1 /im, and 
the critical field Be typically corresponds to $rf < $0. In 
effect, the zero- field conductance minimum is expected to 
be significantly deeper than the other minima, for which 
both electrodes are driven into the normal state. 

In conclusion, we have identified the new transport 
phenomenon in undoped graphene, which manifests it- 
self by periodic magnetoconductance oscillations for the 
Corbino geometry. The relative field-induced conduc- 



tance change reaches experimentally accessible magni- 
tudes AG/Gq > 0.1 for moderate radii ratios. At weak 
doping, the oscillations remain observable for a finite 
range of applied fields [3S]. Additionally, we have pre- 
sented the complete phase diagram in a field-doping pa- 
rameter plane, illustrating the crossover from the field- 
suppressed to the ballistic transport regime, as well as 
the resonances through Landau levels, at which the os- 
cillatory behavior is restored. 

We hope our analysis shall rise some interest in 
Corbino measurements within the graphene community. 
Although the discussion is limited to the system with a 
perfect circular symmetry and the uniform field, partic- 
ular features of the results, including (i) the conductance 
dependce on the total flux piercing the sample area and 
(ii) the formal analogy between dimensionless flux ^d/^o 
and the boundary conditions at zero field, suggest mag- 
netoconductance oscillations should appear in more gen- 
eral situation as well. The work primarily focuses on 
graphene, but the recent study on effective Dirac fermion 
model for HgTe/CdTe quantum wells [26; suggests that 
our findings may also be relevant to such systems. 
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cussions. The support from the Alexander von Humboldt 
Stiftung-Foundation and the Polish Ministry of Science 
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Note added in proof. Recently, I became aware of 
a work on zero-doping Corbino magnetoconductance in 
graphene employing the conformal mapping technique, 
Ref. m\. 
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